Supersymmetric twisting of carbon nanotubes 
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We construct exactly solvable models of twisted carbon nanotubes via supersymmetry, by applying 
the matrix Darboux transformation. We derive the Green's function for these systems and compute 
the local density of states. Explicit examples of twisted carbon nanotubes are produced, where the 
back-scattering is suppressed and bound states are present. We find that the local density of states 
decreases in the regions where the bound states are localized. Dependence of bound-state energies 
on the asymptotic twist of the nanotubes is determined. We also show that each of the constructed 
unextended first order matrix systems possesses a proper nonlinear hidden supersymmetric structure 
with a nontrivial grading operator. 
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INTRODUCTION 

Importance of solvable models in physics is enormous. 
We can acquire qualitative understanding of the com- 
plicated realistic systems by analyzing simplified models 
that grab the essence of a physical reality. These models 
can serve as a test field for approximative methods, or 
can be used as initial solvable systems in a perturbative 
treatment. In this paper, we will focus on the construc- 
tion and analysis of exactly solvable models described by 
the (1 -|- l)-dimensional Dirac equation. 

Such systems lie in the overlap of the quantum field 
theory with the condensed matter physics. The one- 
dimensional Dirac Hamiltonian appears in the study of 
the gap equation of the 1-1-1 dimensional version of the 
Nambu-Jona-Lasinio (chiral Gross- Neveu) model [ij, Q], 
3|, or in the study of the fractionally charged solitons 
4 1 , [H] . It is used in the effective description of the non- 
relativistic fermions: in Q, the Hamiltonian describes 
fermions coupled to solitons in the continuum model of 
a linear molecule of polyacetylene. It is employed in 
the analysis of the quasi-particle bound states associated 
with the planar solitons in superfluid ^He It appears 
in the description of inhomogeneous superconductors [1] 
and in the analysis of the vortex in the extreme type- 
II superconductors in the mean field approximation [9|. 
Last but not least, it is used in the description of carbon 
nanotubes. In the low energy regime, the band structure 
obtained by tight-binding approach can be approximated 
very well with the use of the one-dimensional Dirac op- 
11 1 . The stationary equation fl^l 



erator IIC 



{ia2dj: + Ai{x)ai)(t) = Xcj) 



(1) 



describes dynamics of the low-energy charge-carriers in 
sing le wall carbon nanotubes in presence of magnetic field 



The Green's function (or its spatial trace called diag- 
onal resolvent or Gorkov Green's function) plays an im- 
portant role in the above mentioned systems. It is used 



in solution of the gap equation fl*] or in the extremal anal- 
ysis of the effective action [2] in quantum field systems. 
It is employed in computation of the free energy of the 
inhomogeneous superconductors T3] . It serves in deriva- 
tion of the local density of states (LDOS), the quantity 
that can be measured in carbon nanostructures by the 
spectral tunneling microscopy fl^, [l^]- The results ob- 
tained in this paper will be primarily discussed in the 
latter context. 

The carbon nanotubes are cylinders of small radius 
rolled up from graphene. They can be classified as either 
metallic or semiconducting, in dependence on their elec- 
tronic properties. When no external potential is present, 
the semi-conducting nanotube has a spectral gap which is 
related to a constant value of the potential, Ai = 7^ 0, 
where Py is the value of the canonical momentum in the 
compactified direction. For Ai = = 0, the nanotube 
is metallic as it has no gap in the spectrum. In this case, 
an infinitesimally small excitation is sufficient to move 
the electrons from valence to conduction band. The ac- 
tual value of Py is related to the orientation of the crystal 
lattice in the nanotube, see e.g., [l3|, 

We suppose that the potential Ai(a;) is smooth on the 
scale of the interatomic distance. Otherwise, it would 
be necessary to work with an extended, 4x4, Hamilto- 
nian that would describe mixing of the states between 
the valleys associated with two inequivalent Dirac points 
figt - [2^ . The matrix degree of freedom of in ([T]) is 
the so-called pseudo-spin and is associated with the two 
triangular sublattices that build up the hexagonal struc- 
ture of the graphene crystal; the wave function with ei- 
ther spin-up or -down is identically zero on one of the 
sublattices. 

The inhomogeneous magnetic field can appear due to 
an external source. Alternatively, it can emerge as a 
consequence of mechanical deformations of the lattice. 
Let us make this point clear. Deformation of the lattice 
is described by the vector d = {dx{x,y),dy{x,y)) which 
represents displacement of the atoms in the crystal. The 
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FIG. 1: The nanotube with the twist corresponding to dy --^ x. 
In the untwisted nanotube, the black line would be straight 
(horizontal). 



tion, we briefly review the construction of solvable models 
based on Darboux transformation with focus on the ap- 
plication in the context of carbon nanotubes. Then the 
formulas for Green's function and LDOS of these mod- 
els are provided. We discuss reflectionless systems and 
present two models of twisted carbon nanotubes. The 
last section is left for the discussion. 



associated strain tensor Sij is defined as 



dydy, 



(2) 

The effective Dirac Hamiltonian which describes dynam- 
ics of quasi-particles in the low-energy regime gets the 
form a2{idx + ^2{x)) + (Ti{py + Ai(a;)) -I- lAo, where we 
fixed the Fermi velocity vp = ^- The gauge fields are 
related to the strain tensor ([2]) in this way: A2(a;) = 
{sxx - Syy), Ai(x) = 2sr^y aud Ao(x) = a^^ + Syy up to 



multiplicative constants, see [IJ], [21[, |22| . 

In this context, the potential Ai(a:) in ^ can be in- 
terpreted as the gauge field generated by the twist per- 
pendicular to the axis of the metallic nanotube. The 
angle of the twist d{x) is related to the displacement 
d = (0,/ /S.i[x)dx) by dy{x) = r-d{x) where r is a ra- 
dius of the nanotube. In this way, the constant potential 
Ai{x) — 13 > Q can be associated with a linear displace- 
ment d = (0, jix) that would be generated by the con- 
stant twist illustrated in Figure 1. It opens a gap in the 
spectrum of the metallic nanotube, however, it does not 
confine charge carriers. Indeed, constant potential can be 
understood as a mass term in the Hamiltonian describing 
the free particle. 

In general, the electromagnetic field causes nontrivial 
scattering of the quasi-particles and can even cause the 
appearance of bound states in the system . It is well 
known that the quasi-particles in metallic nanotubes are 
not backscattered by electrostatic potential. This is un- 
derstood as a manifestation of the Klein tunneling [24 1 
and it has been discussed extensively in the literature 
[20| . [25| . It was found recently that the phenomenon 
can be attributed to the peculiar supersymmetric struc- 
ture that relates the Hamiltonian of the system to that 
of the free Dirac particle [l8| . 

Here, we will construct exactly solvable models de- 
scribed by ([T]) where, despite the presence of the effective 
magnetic field, the scattering will be refiectionless and 
the bound states will be confined in the regions where 
the twist gets altered. In the construction, the techniques 
known in the supersymmetric quantum mechanics will be 
employed. We will focus to the spectral properties and 
Green's function of the new systems. The latter one will 
be used for computation of the LDOS. We will provide 
an analytical formula for bound state energies in depen- 
dence on the twist of the nanotubes. 

The work is organized as follows. In the next sec- 



SPECTRAL DESIGN VIA DARBOUX 
TRANSFORMATIONS 

We summarize here the main points of the construction 
of new solvable models which is based on the intertwining 
relations. This scheme is well known in the context of su- 
persymmetric (SUSY) quantum mechanics (26| . There, 
the intertwined second order Schrodinger operators give 
rise to the supersymmetric Hamiltonian while the inter- 
twining operator, identified as the Crum-Darboux trans- 
formation, is associated with the supercharges of the sys- 
tem. In the current case, we will discuss briefiy the tech- 
nique in the context of the first order, one-dimensional 
Dirac equation. We refer to [23| for more details. 

Let us have a physical system described by a solvable 
hermitian Hamiltonian 



h — ia2dx + A 



(3) 



with real and symmetric matrix potential A = A(a;) 
and X extending to the whole real axis. The physical 
eigenstates (solutions complying with prescribed bound- 
ary conditions) form a basis of the Hilbert space. Besides, 
the (formal) solutions of the stationary equation hu = Xu 
are supposed to be known for any complex A. 
We define the operator L by 



L = U-^U-^ = Idx - U'U-\ 
OX 



(4) 



where U' — dU/dx. The matrix U — (ui,M2) is a chosen 
solution of the equation hU — U A where the matrix 
A — diag{Xi, \2) has fixed real elements. The vectors 
Mi(2) satisfy /iui(2) = Ai(2)Wi(2) and are chosen to be 
real. They form the kernel of L, LU = and do not need 
to be physical. Next, we define the hermitian operator h 
with the potential term explicitly dependent on ui and 
1*2 and corresponding eigenvalues Ai and A2, 



h = 



h + i[a2, U'U-^] = a2hn2 + (J2[cj2, C/AC/-1] 

/ uJaiU2 ula^U2 \ , . 



We used here the identity 

Idx = -i(T2{h - A), 



(5) 



(6) 



which will be employed extensively in the following text. 
Notice that as long as ui and U2 correspond to the same 
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eigenvalue Ai = A2, h reduces to a unitary transformed 
seed Hamiltonian, h = a2ha2, for any A. The Hamiltoni- 
ans ([3]) and ([5]) satisfy the fohowing intertwining relations 
mediated by L and L\ 

Lh = hL, L'^h = hLK (7) 

The conjugate operator can be written as = —dx + 
V'V~^ where V — (U^)^^ = (vi,V2)- The columns vi 
and V2 satisfy i^Ui(2) ~ and are solutions of hvi{2) — 
\(2)Vi{2)- There holds 

hV = V K. 

Each of the equations {h~-\)ip — and {h—X)(p = has 
two independent formal solutions, let us denote them ipx^ 

and respectively. For A 7^ ^1(2), the operators 

L and work as one-to-one mappings between the two 
subspaces spanned by ^\ and ip\, They transform 
the (formal) eigenvectors of h into the formal eigenvectors 
of h and vice versa. 

Let us consider now the four-dimensional subspace 
spanned by the solutions of {h — Ai(2))</? = 0. Two of 
the solutions, the vectors ui and U2, compose the ma- 
trix U. We can use the other two vectors to define the 
matrix U_ = {Hi, 112)1 which satisfies hU_ — U_A but is 
not annihilated by L. Similarly, we can define the ma- 
trix y_ = (w]^, W2) from the solutions of {h — Ai(2))'^ = 
which satisfies hy_ = V.A, but no linear combination of 
Vi and V2 is annihilated by L'^ . The intertwining opera- 
tors then transform the matrices as 
Hence, we get L'<LU = L'LU ^ Wv = WV_ = 0. 
The latter equalities can be understood as the implica- 
tion of the alternative presentation for the products of 
the intertwining operators, 

LL^ = 01- \i){h - A2), L^L ={h- Xi){h - A2). (8) 

The spectrum of h is identical with the spectrum of h 
up to a possible difference in the energy levels Ai and/or 
A2. These energies are in the spectrum of either h or 
h if and only if the associated eigenvectors comply with 
the boundary conditions of the corresponding stationary 
equation. We will discuss specific examples where the 
spectrum of the new Hamiltonian h contains additional 
discrete energies that are absent in the spectrum of h. 
The eigenvector 0^ of h corresponding to the energy level 
7^ Ai, A2 can be expressed in terms of L and the 
eigenvectors 4>k of h (jKpk = Xk4>k), 

^^= V(A.-Ar)k-A2) ' 

When defined in this way, the probability densities of (j)k 
and 0fc coincide. 

It is worth noticing that the system h inherits integrals 
of motion of h. Indeed, if S commutes with h, then the 



operator S = LSL^ generates a symmetry of h, [h, S] = 
0. We will discuss this point in more detail in the context 
of the reflectionless models. 

The potential term of /i in ([5]) ceases to have a direct 
interpretation in the context of carbon nanotubes with 
the radial twist. As we are interested in the analysis of 
namely such systems, we require h to be equivalent to 
the Hamiltonian in ([1]); the term proportional to either 
(7i or (T3 in (lU should vanish. As these coefficients depend 
both on the potential of the seed Hamiltonian h and on its 
eigenvectors, it is rather difficult to meet this requirement 
in general. Instead, let us consider two special cases. 

First, let us fix the initial Hamiltonian as 

hi = i(T2dx + m(T3 + Aicri, (9) 

where m > 0. We take Ai = m and A2 = and de- 
note Ui = U = (ui,U2) and Vj = V = (1^1,^2), where 
explicitly ui = (uii,0)^, U2 = (ui2,U22)^ and 

\ / \ U11U22 «22 / 

Comparison of the two matrices tells that if ui (or U2) 
is a bound state of h, then vi (or V2) cannot be bound 
state of h. Vice versa, if vi (or V2) is a bound state of 
h, then ui (or U2) is not normalizable. Using ([5|) and 
uJaiU2/ det [// = 1, we get the Hamiltonian hi, 

/i/ = icr29j; - ( Ai + m— ) (71, (11) 
V "22/ 

with the required form of the potential. 

In the second case, we take the seed Hamiltonian as 

hii = i(j2dx + (Ai + 771)0-1 (12) 

and fix Ai = — A2 > 0. The vectors Ui(2) are chosen as 
ui = (uii,7i2i)"^ and 77,2 = cratii. They satisfy u[aiU2 = 
0. The matrices Un = U and Vn = V are in this case 

Uii = ( V Vii = \( ""1 ) (13) 

\^ 7*21 -7i21 / 2 \ 7t2i -7*21 / 

and the Hamiltonian ([S]) acquires the form 

hii = ia2dx~ fAi-|-m-A /"^'''^ Vi- (14) 

We can deduce that if 7ti is a bound state of h, so is the 
vector U2 and neither 7;i or V2 can be normalized. Vice 
versa, if vi and 772 are bound states of /i, the vectors ui 
and U2 are not normalizable. 

In the end of the section, let us notice that there is an 
alternative interpretation in dealing with the intertwining 
relations and the involved operators. Inspired by the 
SUSY quantum mechanics, we can define the extended, 
first order matrix operators 

«^(S2). 0)^ e-(-°7.0 

(15) 
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which establish the N = 2 (nonhnear) supersymmetry. 
The grading operator F = diag(l,— 1) classifies the 
Hamiltonian T-L as bosonic ([H,r] ~ 0), while both Qi 
and Q2 = «rQi are fermionic, {Qa, F} = for a = 1, 2. 

Contrary to the SUSY quantum mechanics based on 
the second order matrix Hamiltonian, here both the 
Hamiltonian T-L and the supercharges Qi^2 are the first 
order differential operators. The associated (nonlinear) 
superalgebra 

[H,Qa]=0, {Qa,Qb} = 25ab{n-\l){n-\2) (16) 

encodes the intertwining relations ([7]) together with the 
factorization (IH). 



GREEN'S FUNCTION AND LDOS FOR THE 
TWISTED NANOTUBES 

We shall derive formula for the Green's function of h 
in terms of the intertwining operator L and the Green's 
function of the initial Hamiltonian h. In the end of the 
section, we will discuss the explicit form of the LDOS for 
the systems described by hi and hn of the form (fTTj) and 
(ITil) corresponding to the twisted carbon nanotubes. 

Let us start with the hermitian Hamiltonian h — 
i<j2dx + A(a;). The potential term is required to be real 
and symmetric. The (generalized) eigenstates of h 
have to satisfy the following boundary conditions 



>±oo 



/±(A,x), Ae 



(17) 



The symbol ^ means here that the elements of the eigen- 
vector are proportional asymptotically to the function 
f±{x, A). We prefer to leave the boundary conditions un- 
specified explicitly at the moment. They will be discussed 
for the refiectionless models later in the text. 

The Green's function associated with the Hamiltonian 
h is defined as a solution of the equation 



{h-X)G{x,y-\)^5{x~v), A G C. 



(18) 



It has to satisfy the same boundary conditions as 
the eigenstates of /i, i.e. the matrix elements of the 
Green's function are proportional to /±(A, x) in the limit 
X — >■ ±oo. Being effectively the inverse of {h — A), the 
Green's function is not well defined for A from the spec- 
trum a{h) of h. It has simple poles for A corresponding to 
discrete energies. If A is in the continuous spectrum, then 
we can find the limit G'^(x, y; A) — lim^_).o G(a;, y; AztiTy), 
see e.g. 

The differential equation in p7|) has two formal inde- 
pendent solutions iIj\{x) and £,x{x) for any A G C. For 
A ^ cr{h), we can fix ip\ and such that each of the 
functions complies with the boundary condition in one 
of the boundaries; i.e. we fix i/'A(a^)|a;->-+oo ^ f+ix^X) 



and ^\{x)\x^-oo ^ f-{x, A). These functions can be em- 
ployed in the construction of the Green's function in the 
following way 



G(x,y; A) = 



Mx)UyV0{x ~y) + U^)Myfo{y - x) 



(19) 



where is the step function. The quantity 

W{^,0=i^{^f<y2^{^) (20) 



is the analog of Wronskian for Dirac equation. It is 
constant for two independent solutions t/'a and cor- 
responding to the eigenvalue A of h. Indeed, direct cal- 
culation with the use of © shows that dxW{i}},C) — 0. 
The Green's function defined in (|19l) then solves (ITS)) and 
manifestly satisfies the prescribed boundary conditions 
for X — >■ ±cxD. 

Let us pass to the system described by h and con- 
struct its Green's function with the use of ([T^ . We 
suppose that L transforms appropriately the boundary 
conditions associated with h to the boundary conditions 
prescribed for the eigenstates of h. We can define the 



functions ^a 



and ^A 



^(A-Ai)(A-A2) ^ ^(A-Ai)(A-A2)' 

They solve hxjjx = AV'a, h^^x = A^a and satisfy the pre- 
scribed boundary condition in -|-oo or — cx), respectively. 
The Green's function associated with h can be written 
then as 



G{x,y\X) 



M^)Uy)^e{x - y) + Ux)i^x[yVe{y - x) 



1 



Wi^x.ix) 



T^(V'A,6) 
{Li;x){x){LixV{y)e{x-y) 
(A- Ai)(A-A2) 



{Lix){x){L^Jxf{y)e{y-x) 



(A- Ai)(A- A2) 



(21) 



We used the fact that the Wronskian is invariant with 
respect to the Darboux transformation (|4]), W{il)XiS,x) = 
W{il)x,£.x)- We refer to [3 or [H] where the proof of this 
relation can be found. Let us mention that the a dif- 
ferent supersymmetric approach to Green's functions of 
Dirac operators was examined in 'sci] where a modifica- 
tion of the standard supersymmetry (based on second- 
order Hamiltonians) was discussed. 

The eigenvectors of h can be written as 



itx^ C{X,x)ipx, C{X,x) 



-1(72- 



X-U{x)AU-'^{x) 



v/(A-Ai)(A-A2) 
(22) 

where we used (O again. This allows us to rewrite the 
Green's function (1211) in particularly simple form 



G{x,y;X) ^ C{X,x)G{x,y;X)£'^{X,y). 



(23) 



Hence, G{x, y; A) can be obtained by purely algebraic 
means without the use of any differential operator; it 
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can be obtained just by multiplication of G{x, y; A) with 
simple matrix operators (|22|) . 

The local density of states p{x, A) associated with h is 
computed in the following manner 



p{x, A) = lim 

TT ImA^04 



linTrG{x,x;X), (24) 



where the trace is taken over the matrix degrees of free- 
dom. Using (1231), we can write LDOS p for h as 



Im Tr {C{X, xfC{\ x) G(x, x; A)) 



p(a;, A) = lim 

ImA^O- 

(25) 

Notice that the formulas ([23l) and ([25]) are valid for a 
general class of the seed Hamiltonians with real and sym- 
metric potential. 

In the literature (see, e.g. [H, 0, [ll]), the operator 
G{x, x; A) is called Gorkov Green's function or diagonal 
resolvent of h. The Green's function of the Schrodinger 
operators and generalized Sturm-Liouville equation was 
studied in 31 1 and (s^ in the context of intertwining 
relations. 

We turn our attention to the systems represented by 
hj and hn which describe the carbon nanotubes with the 
radial twist. It is supposed that the Green's functions of 
both hi and hu are known. We denote them Gi{x, y; A) 
and Gii{x,y;\). The operators Li{\,x) and Lii(\,x) 
based on Ui and Un respectively acquire particularly 
simple form 



and 



£/(A,x) 



Cii{\,x) 



1 







-A 



-m + A 

M22 



Ai^^ -A 

J- till 

A -Ai2^ 

U21 



The trace of the G/ can be computed directly in terms 
of the vectors ui and U2- A straightforward computation 
gives 



Tr(G/(a:,a:;A)) = 50 



2 t 
m u[ui 

2A(A-m)(detC//)2 



m m2 5o(ul'«i)(u2U2) 

-ff3 



A — to' 



2A(A- m)(detC//)2 



A — TO i 

51 ul^aiU2 



(26) 

Here we used the abbreviated notation go = 
Tr G/(x, x; A) and = Tr(crjG/(a;, x; A)) for j = 1,3. 
We can obtain similar expression for the trace of the 
G//(a;,a;; A): 



Tr(G//(x,a:; A)) = go + 



2\lgo{u\uif 
(A2-A2)(detC///)2 



2\\u\ui 



(A2-A2)(detC/„)2 



t A t ^ 

- gi—u[(Tiui 



(27) 

The notation used here is like in ([26]) with the replace- 
ment of Gi{x, x; A) by Gii{x, x; A). 



PERFECT TUNNELING IN THE TWISTED 
CARBON NANOTUBES 

There exists an exceptional class of exactly solvable 
systems whose Hamiltonian h is intertwined with the 
Hamiltonian of the free particle. The peculiar and sim- 
ple properties of the latter model are manifested in these 
systems as well. In particular, they share the trivial scat- 
tering characteristics of the interaction-free model, i.e. 
they are reflectionless. The eigenstates of both the free- 
particle system and the reflectionless models are subject 
to the same boundary conditions; the scattering states 
have to be oscillating in the infinity while the bound 
states should decay exponentially for |x| — )■ 00. Addi- 
tionally, the reflectionless systems inherit the integral of 
motion that in the free particle system plays the role of 
generator of translations. 

The stationary equation hcj) — Xcj), where h = ia2dx + 
mas, is translationally invariant, i.e. the Hamiltonian 
commutes with p — ~idx- We can find the common 
eigenstates of h and p. The latter operator distinguishes 
the two scattering states corresponding to each doubly 
degenerate energy level. It annihilates the singlet states 
u+ = (1,0)"^ and u_ = (0,1)^ that correspond to the 
edges A = ±m of the positive and negative part of the 
continuous spectrum (which are called the conduction 
and the valence band respectively in the context of nan- 
otubes). The involved operators close the nonlinear su- 
peralgebra 



[p,h]^0, {p,p} = 2ih-m){h + m), 



(28) 



which is graded by the parity operator F = Ra^ {RxR = 
—X, = 1). Let us stress that this supersymmetric 
structure is completely different from (|16p . In this case, 
the supersymmetry is rather hidden; the two fold de- 
generacy of energy levels, distinguished by the integral 
of motion p, emerges within the spectrum of the unex- 
tended Hamiltonian h. 

The Hamiltonian h inherits a modified version of the 
nontrivial integral of motion p. It can be found by dress- 
ing of the initial symmetry operator. 



P 



LpL\ [p,h]^Q. 



(29) 



It annihilates the states vi and V2 together with the vec- 
tors v± which are defined as v± = Lu± . The operator p, 
like p in the free particle model, reflects the degeneracy of 
the spectrum; it can distinguish the scattering states cor- 
responding to the same energy level. The superalgebra 



6 



(l28)l can be recovered in the modified form 

[p,h]=0, {p,p}^2Ch^^m^)Ch-X,fCh-X2f. (30) 

Hence, tlie square of p is the spectral polynomial of h. 
It is worth noticing that the same algebraic structure, 
the hidden supersymmetry, was discussed in detail for 
both relativistic and nonrelativistic finite-gap systems in 
[ssf . [sl j. [ssf . [3^. In this context, the integral p can be 
identified as the Lax operator of the system represented 
hyh. 



Single-kink system 

The first model will be derived with the use of the seed 
Hamiltonian hj in (|9]) with Ai(a;) — 0, 

hj = iCT2<9x + mas. 

We will compute its LDOS and discuss the realization of 
the parity operator of the hidden supersymmetry. 

We require that the new Hamiltonian h has a single 
bound state with zero energy. To meet this requirement, 
we fix the matrix Ui as 



2 / 1 — sinh 1 



\ m \ coshmx 
and the intertwining operator L as 

Li = -i-d, + m(^^Q -tanhmx 
The explicit form of the matrix Vj is then 



(31) 



1 







V = - 

2 V tanh mx sechma; 



The associated Hamiltonian hi then reads 

hi = i(J2dx + mai tanhma;. (32) 
The operator hi has the normalized bound state V2 



FIG. 2: The metallic nanotube with the twist associated with 



nanotube, the black line would be straight. 



The presence of the step function 6 reflects that fact that 
imaginary part of (I26|) for |A| < m is zero and, hence, 
p{x, A) vanishes identically. The formula (I33p can be 
rewritten with the use of the LDOS of the free particle 



Piix,X) 



|A| 



:^(A2 -m2) 



and the density of probability of the bound state V2 , 



771 4- 

pi{x, A) = pi{x, A) ( 1 - J^2'"2'^'2 



(34) 



The coefficient of the second term is just the difference 
of the densities of states of h and h. 



{pi - pi)dx 



to6I(A2 -to2) 



7r|A|v/|A2 -to2 



Let us notice that the difference of densities of states 
for Dirac particle on the finite interval with Dirichlet 
boundary conditions was discussed in i29] . 

The hidden superalgebra (PH)) . closed by hi and pi = 
LipL\, reads explicitly 

[hi,pi] ^ 0, {pi,pi} = 2{hi - mfh]{hi + m). 

The parity operator F — Ra^, = 1, classifies hi 
and Pi as, respectively, bosonic and fermionic operators, 
[hi,f] = {pi,f} = 0. 

The potential in (|32l) can be associated with the dis- 
placement vector d — (0, In cosh mx). The corresponding 
twist of the metallic nanotube is illustrated in Figure [2] 
Hence, the nanotube is twisted in one direction up to 
the center (origin) where the orientation of the twist gets 
changed. 



^2 — { 0, \l sech mx 



Let us notice that the operator (|32|) appears in descrip- 
tion of many physical systems, e.g. in the continuum 
model for solitons in polyacetylene @ or in the analysis 
of the static fermionic bags of the Gross-Neveu model Q . 

We can use together with to compute the 
LDOS of the system. It acquires the following simple 
form 



piix,X) 



2|Ap — m^sech^ma: 
27r|A|y|A2^^^ 



e{X'^-m^). (33) 



Double-kink model 

Here we construct the system with two bound states. 
We shall employ the scheme discussed in (IT^ - p^ . Fix- 
ing Ai(a::) = in (fT2|). we get the free particle Hamilto- 
nian 

hii = i<J2dx + mtTi. 
We choose the components of Un as 

uii — —j= coshfcx, U21 = —j= cosh(fca; -I- a) 
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where 



1 , m — k 
a = - log r, k 



m? — \\, < Ai < m. 



2 °m + k 

The intertwining operator acquires a diagonal form 

(35, 

The formula then provides the explicit form of the 
Hamiltonian hn 

~ / , cosh^ kx + cosh^ffcx + a) ' 

hii ~ lO-20x + ^JTl + Ai — ■ ^ I (Ti 



cosh kx cosh(fca; + a) 
= i(T2dx + (ti — fctanhfca; + fctanh(fca:; + a))(Ti.(36) 

The potential term is asymptotically equal to mai. The 
system has two bound states represented by the normal- 
ized vectors vi and V2 = vi where 

Vk 

vi ~ -^(sechfca;, sech(fca; + a))'^ . 



Notice that the equation (l36l) appeared in the analy- 
sis of the Dashen, Hasslacher, and Neveu kink-antikink 
baryons in Gross-Neveu model js^l • 

The local density of states in the current system can 
be computed directly with the use of (P7|) . We get 



Pii{x,X) = 



|A| ^1 — 2(\^-X^) (sech^fcx + sech^(fcx + a)) 



7r^|m2 - A2| 



(37) 



Likewise in the preceding example, it can be written as 
the LDOS of the free particle corrected by the term pro- 
portional to the probability density of the bound states, 

~Ph{x,\) ^ Pu{x,\) (^^- (^^'^2) ^ ' (38) 

where pii{x,X) = pi{x,X). This time, the difference of 
the densities of states is 



ADOS' = [ {pq- pi)dx = - 



2fc|A|6i(A2 - m^) 



rv/|m2- A2|(A2 - A?) 

The hidden superalgebra associated with the system, 

[hii,Pii] = 0, {pii,Pii} = nhjj - m^)ihjj - Xlf, 

is graded by the operator T = RRaUi where Raf{x) = 
fix + a)Ra, RaR = RR^a i?(-Ra)"^ and a = -|. 
This grading operator (represented in another form) was 
also discussed in [sl]. 

The vector potential in ([BG)) corresponds to the dis- 
placement dy = mx — Incoshfccc -|- lncosh(fcx -f a). The 
corresponding twist of the metallic nanotube does not 



FIG. 3: The metallic nanotube with the twist associated with 
dy ~ mx + In ^jjjj ^Yie Hamiltonian 



FIG. 4: The nanotube associated with the Hamiltonian ([36]) 
and the twist corresponding to (|39p . The constant part of the 
magnetic field in p6|) . Aa/g = rn, can be attributed to the 
external magnetic field or to the semi-conducting character of 
the nanotube. 



change its orientation asymptotically, see Figure [3] for 
illustration. 

We can find another physically interesting setting de- 
scribed by hjj. We can divide the vector potential 
in (j36p into two parts. The first part is associated 
with the asymptotically vanishing twist of the nanotube. 
At = —ktanhkx + fctanh(fca:; -I- a). The second part is 
constant, Amg = Jti, and corresponds to the homoge- 
neous external magnetic field which is parallel with the 
axis of the nanotube. Hence, hjj describes the metallic 
nanotube which is asymptotically free of twists, however, 
the external constant magnetic field is present. See Fig- 
ure [5] for illustration. 

The uniform external field opens a gap of the width 2m 
in the spectrum while the asymptotically vanishing twist 
induces two bound states in the gap. The model allows 
to compute the bound state energies as a function of an 
asymptotic (global) twist. Indeed, the twist associated 
with At is 



dy = In 



cosh(fca; -I- a) 
cosh kx 

The asymptotic twist corresponds to 



(39) 



Sd — I lim dy— lim dy \ = 2\a 

x^oo a;— > — oo 



In- 



— m? — \\ 

m + y/ m? — \\ 
(40) 

The dependence of the bound state energies on 5d then 
acquires the following simple form 



Ai = ±2m- 



e 2 



(41) 



and is plotted in Figure [S] 

Up to now, the twisted nanotubes were considered 
to be metallic. The analysis can be extended to semi- 
conducting nanotubes without any difficulties; a con- 
stant, nonzero, part of the vector potential Ai has to 
be associated with the internal characteristics (the orien- 
tation of the hexagonal lattice) of the nanotube. Let us 
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Sd 



FIG. 5: The spectrum of the Hamiltonian ((36]). The asymp- 
totic twist of the nanotube (|40p induces bound states of en- 
ergies (|41[) . The parameter m is proportional to the inverse 
of the radius of the nanotube, see |43l ]. 



notice in this context that the metalhc nanotube can be 
converted into the semi-conducting one just by switching 
on the constant magnetic flux parallel to the axis of the 
nanotube. This fact was experimentally confirmed in (38| 
and coined as Aharonov-Bohm oscillations of the carbon 
nanotubes. Turning back to (|36p . we can interpret the 
Hamiltonian as the energy operator of the semiconduct- 
ing nanotube with a radial twist associated to At- The 
constant part py = Amg of the potential appears due to 
the semiconducting nature of the nanotube. 



DISCUSSION 



The expressions (|34p and 
unified form 



can be written in the 



1-E 



-^1 



(A2 



where the sum is taken over the normalized bound states 
of hj^jj) annihilated by Lj(^jjy It manifests a decrease 
of the LDOS in the regions where the bound states are 
localized. However, it is rather just a peculiar property of 
the discussed reflectionless models [39|. In general case, 
the LDOS ([27l) of hji acquires the following form in terms 
of the vectors vi and V2, 

TT{Gn{x,x;X)) = Tt{Gii(x,x] X)) 



2Xlvlvi (^govlvi + gsvlasvi 



When hjj is equal to the free particle Hamiltonian, the 
coefficient of v\vi reduces to a constant. Nevertheless, 
this apparently does not hold in the general case. 

We restricted our consideration just to the systems de- 
scribed by the Hamiltonian h = ia2dx + Ai{x)cri. How- 
ever, the potential term of the seed Hamiltonian in ([3]) 
can acquire quite generic form, yet keeping valid the for- 
mulas (l23l) and (|25|) for the Green's functions and for the 
LDOS. Other results are more sensitive to the explicit 
form of the potential. The term A2{x)a2 cannot cause 
any substantial modifications; it would play just the role 
of non-physical gauge field. In contrary, impact of the di- 
agonal term IAq -I- A3{x)a3 in h would be much deeper: 
in general, it would break the symmetry a{h) = —a{h) of 
the spectrum a{h) of h. In the context of Dirac particles 
in the carbon nanotubes, the potential A3(x)cr3 would 
have different sign for the spin -up and -down compo- 
nents of wave function, i.e. this potential would change 
the sign on the two sublattices that form the crystal. 
Physical realization of such a scenario in the considered 
condensed matter system is not clear. It is remarkable 
that the Darboux transformation (U) does not alter the 
form of Aq ; the new Hamiltonian h shares the same elec- 
trostatic potential as the seed Hamiltonian. We notice 
that the electrostatic potential can be also altered via the 
so-called 0-th order supersymmetry, as it was discussed 
in 



18| 



(A2-A^)(det Vy 



In the discussed systems represented by the station- 
ary equation ([1]), the analysis of the bound states can be 
facilitated by the fact that the square of the Dirac Hamil- 
tonian takes the form —d^ + A\ + cr^A'i. The existing 
tools (see, e.g. [i^) for the analysis of the Schrodinger 
operators can be exploited to reveal spectral properties 
of the Dirac Hamiltonian. In this context, let us men- 
tion that interesting results were obtained by the spectral 
analysis of general class of deformed quantum waveguides 
described by Schrodinger equation 4l|. We believe that 
similar analysis for the carbon nanostructures described 
by the one- or two-dimensional Dirac Hamiltonian would 
be fruitful. 

The presented analysis is qualitative. The equation 
dl]) is a good approximation for the quasiparticles in car- 
bon nanotubes only for small region of the momentum 
space where the linear dispersion relation is valid. When 
the gap opened by the pseudo-magnetic field in the spec- 
trum is too big, nonlinear (the so-called trigonal warping) 
terms have to be included into the Hamiltonian. In 
the article, we neglected surface curvature of the nan- 
otubes. The tubular surface prevents the 7r-orbitals of 
the carbon atoms to be parallel to each other. This im- 
plies presence of additional pseudo-magnetic fields in the 
Hamiltonian. However, in case of armchair nanotubes, 
these additional gauge fields can be transformed out Q . 

The examples presented in the text suggest that the 
non-uniform radial twist can induce bound states in the 
nanotube. In particular, the second model with double- 
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kink potential provides an interesting qualitative insight 
into realistic experimental setting: the nanotube with 
asymptotically vanishing twist is immersed into the ho- 
mogeneous magnetic field. Besides the explicit formula 
([55)1 for LDOS, the model predicts the appearance of 
bound states and the formula (|^T|) fixes their energies 
in dependence on the asymptotic twist. The model can 
be simply tuned with the use of perturbation techniques. 

The supersymmetry can be very useful for construc- 
tion of the models with more complicated (yet asymptot- 
ically constant) twist inducing richer spectral properties. 
The formalism presented in the second section can be re- 
peated to produce a chain of solvable Hamiltonians, h, h, 

h, by taking the last constructed operator as the seed 
Hamiltonian for the new system. These new solvable sys- 
tems shall provide insight into the deformation-induced 
spectral engineering of carbon nanotubes. The reflection- 
less models are particularly important in this context; 
they are analytically feasible and possess nontrivial (su- 
per)symmetry, analog of (|29|) and (|30|) . The considered 
double-kink example suggests that the number of bound 
states could be in a simple relation to the vector potential 
of the Hamiltonian; the number of bound states might be 
proportional to the number of minima of the potential. 
Verification of this hypothesis goes beyond the scope of 
the current paper and should be discussed elsewhere. 
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